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Abstract. Forced turbulence simulations are used to determine the turbulent kinematic viscosity, v%, from the 
decay rate of a large scale velocity field. Likewise, the turbulent magnetic diffusivity, r^t, is determined from the 
decay of a large scale magnetic field. In the kinematic regime, when the field is weak, the turbulent magnetic 
Prandtl number, Vt/rjt, is about unity. When the field is nonhelical, r/t is quenched when magnetic and kinetic 
energies become comparable. For helical fields the quenching is stronger and can be described by a dynamical 
quenching formula. 
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1. Introduction 

The concept of turbulent diffusion is often invoked when 
modeling large scale flows and magnetic fields in a tur- 
bulent medium. Turbulent magnetic diffusion is similar to 
turbulent thermal diffusion which characterizes the turbu- 
lent exchange of patches of warm and cold gas. This con- 
cept is also applied to turbulent magnetic diffusion which 
describes the turbulent exchange of patches of magnetic 
field with different strengths and direction. Reconnection 
of magnetic field lines is not explicitly required, but in the 
long run unavoidable if the magnetic power spectrum is to 
decrease toward small scales. The idea of Prandtl is that 
only the energy carrying eddies contribute to the mixing 
of large scale distributions of velocity and magnetic field 
structures. This leads to a turbulent magnetic diffusion 
coefficient r)t ~ \Ui, where U is the typical velocity and I 
the scale of the energy carrying eddies. For the kinematic 
turbulent viscosity one expects similar values. Analytic 
theory based on the quasilinear approximation also pro- 
duces similar (but not identical) values of r) t and v t (e.g. 
Kitchatinov et al. 1994). 

It is usually assumed that the values of rjt and u t are 
independent of the molecular (microscopic) viscosity and 
magnetic diffusivity, v and r\. However, in the context of 
the geodynamo or in laboratory liquid metals the micro- 
scopic magnetic Prandtl number, P m = v/rj is very small 
(« 10~ 5 ). This has dramatic consequences for the mag- 



netorotational instability (Balbus & Hawley 1991). This 
instability is generally accepted as the main mechanism 
producing turbulence in accretion discs (Balbus & Hawley 
1998). For sufficiently small values of P m , however, this 
instability is suppressed (Riidiger & Shalybkov 2002). On 
the other hand, the Reynolds number of the flow is quite 
large (10 5 . . . 10 6 ) and the flow therefore most certainly 
turbulent. This led Noguchi et al. (2002) to invoke a tur- 
bulent kinematic viscosity, v t , but to retain the micro- 
scopic value of rj. The resulting effective magnetic Prandtl 
number they used was 10 -2 - big enough for the mag- 
netorotational instability to develop. On may wonder, of 
course, why one should not instead use turbulent values 
for both coefficients, i.e. v t /rjt w 1. This would lead to 
even more favorable conditions for the magnetorotational 
instability (Riidiger et al. 2002). 

Similar constraints have also been reported for the 
convection-driven geodynamo: Christcnscn et al. (1999) 
found that there is a minimum value of P m of about 0.25 
below which dynamo action does not occur at all. Similar 
results have also been reported by Cattaneo (2003). These 
results are disturbing, because both for the sun and for 
the earth, P m <C 1. For P m of order unity, on the other 
hand, earth-like magnetic configurations can more easily 
be reproduced (see Kutzner & Christensen 2002). 

Because of these restrictions, one wonders whether the 
effective magnetic Prandtl number to be used is not P m , 
but rather the turbulent value, P m ,t = ^t/vt- This raises 
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the important questions whether P mit is actually of order 
unity and whether it is independent of the microscopic 
value, P m . The aim of this paper is to estimate the value 
°f Pn,t using turbulence simulations. 

The knowledge of the value of P m> t is also important 
for the solar dynamo. The qualitative properties of the dy- 
namo depend on the relative importance of the large scale 
flows and hence on the magnitude of rjt ■ If r) t is too large, 
the influence of a meridional flow of say 10 m/s is small so 
that only little modification can be expected for the basic 
afi-dynamo (Roberts & Stix 1972). In this case, however, 
we know that conventional dynamo models of the solar 
activity cycle have difficulty to explain Sporer's law of 
equatorward sunspot migration. The alternative that the 
resulting poleward migration can be overcompensated by 
an internal equatorward flow requires a sufficiently small 
value of 7ft, which implies that P mjt > 1 (Choudhuri et al. 
1995; Dikpati & Charbonneau 1999; Bonanno et al. 2002). 

Given the importance of the value of the turbulent 
magnetic Prandtl number it is useful to assess the problem 
using three-dimensional simulations of turbulent flows. We 
determine v t and rj t by measuring the decay rate of a 
large scale (mean) velocity and magnetic field, u and B, 
respectively. We emphasize that we are not addressing the 
question whether vt and rjt can really be used in studies 
of the dynamo or the magnetorotational instability, for 
example. 

We consider weakly compressible nonhelically forced 
turbulence and use a model similar to that of Brandenburg 
(2001), but with kinetic helicity fluctuating about zero. 
Dynamo action for such a model has recently been con- 
sidered by Haugen et al. (2003), but it sets in only at mag- 
netic Reynolds numbers above ~ 30, which is not the case 
in the present simulations. We begin however by first re- 
viewing the basic results for the values of vt and rjt within 
the framework of the quasilincar (Roberts & Soward 1975, 
Riidiger 1989) and other approximations. 

2. Results from quasilinear approximation 

For steady homogeneous isotropic turbulence the correla- 
tion tensor is independent of x and t, i.e. 



(uj(x, t)u'Ax + i, t + r )> = Quit, t), 



(1) 



where angular brackets denote an ensemble average and 
primes fluctuations about the average. In the quasilinear 
approximation the transport coefficients are conveniently 
expressed in terms of the Fourier transformed correlation 
tensor, Qij(k,u), which is normalized such that 



(2) 



For the turbulent viscosity and the turbulent magnetic 
diffusivity one finds respectively (Riidiger 1989) 
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v 6 k b Qu(k,Lu) 
(w 2 + ^fc 4 ) 2 



dk cLo, 



(3) 



•ft = 3 



rjk 2 Q H (k,m) 

LU 2 + T] 2 k 4 



dk dio. 



(4) 



Obviously, both quantities are of the same order of mag- 
nitude, but they are not identical. In the limits ^, — > 
the expressions are drastically simplified, i.e. 



v% = 
and 

Vt = 



t_xj 

J (u'(x,t) -u'(x,t + T))dT 

— oo 

OO 

~ J (u'(x,t)- U '(x,t + T))dT, 



(5) 



(6) 



0.4. 



(7) 



so that for the turbulent magnetic Prandtl number is 
v t 2 

Pn,t — — — — 

?7t 5 

This results is similar to that of Nakano et al. (1979) for 
the thermal Prandtl number. 

Riidiger (1989) lists a number of other approaches for 
calculating turbulent transport coefficients, which all yield 
Prandtl numbers around or below unity. One particular 
approach is the renormalization group analysis which was 
applied to turbulence by Forster et al. (1977) for the case 
of a passive scalar, and later by Fournier et al. (1982) to 
the case with magnetic fields. These results are valid in the 
long-time large-scale limit, and the value of P mj t turned 
out to be close to 0.7; see Eq. (23) of Fournier et al. (1982). 

Kitchatinov et al. (1994) use a mixing length approxi- 
mation where terms of the form d/di— z^V 2 are replaced by 
t~1., where r rnrr is the correlation time of the turbulence. 



They find v t = (4/15)r, 



corr"*rms 



and rj t = (l/3)r, 



corr »rms 1 



so P mj t = 4/5 = 0.8. Yet another approach is the r- 
approximation where triple correlations are replaced by 
a damping term that is proportional to the quadratic 
moments (e.g. Kleeorin et al. 1996, Blackman & Field 
2002). Here no Fourier transformation in time is used. 
This gives, as before, rjt = (l/3)ru 2 ms (where r is now in- 
terpreted as a relaxation time), but v t = (2/15)rit 2 ms , so 
Pm,t = 2/5 = 0.4. This is half the value obtained from the 
mixing length approximation, but the same as in Eq. (J7J- 

The fact that in all these cases P mt is less than unity 
can be traced back to the presence of the pressure term 
in the momentum equation. If this term is ignored (as in 
pressureless Burgers turbulence or 'burgulence') one al- 
ways gets P m .t = 1. 

It is tempting to speculate that the discrepancy be- 
tween the different analytic approaches is related to the 
validity of some idealizing assumptions made in order to 
apply the quasilinear and other approximations. Clearly, 
additional approaches are needed to get a more complete 
picture regarding the correct value of P mi t- It is neverthe- 
less encouraging that P m! t does not strongly deviate from 
unity. 

In the remainder of this paper we estimate vt and rjt 
numerically by considering the decay of an initial large 
scale velocity or magnetic field, respectively, in the pres- 
ence of small scale turbulence. 
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3. The model 



The equations describing compressible isothermal hydro- 
magnetic flows with constant sound speed, c s , are 



Dm 9 _ , J x B 

_ = -c s 2 Vmp + -— + J 

Dlnp _ 
OB 

— = V x (u x B) + nV 2 B, 
at 



(8) 



(9) 



(10) 



where u is the velocity, p the density, B is the magnetic 
field, and J = V x B/^q is the current density with /io 
being the vacuum permeability. The viscous force is 



3 v v 



u 



2S • Vlnp) , 



(11) 



|<%V • u is the traceless rate 



where S v — ^y^,] 
of strain tensor 

We solve the equations using the Pencil Code 1 , which 
is a memory-efficient sixth-order finite difference code us- 
ing the 2iV-RK3 scheme of Williamson (1980). For most 
of the simulations a resolution of 128 3 meshpoints is used, 
but in Sect. 0a higher resolution of up to 512 3 meshpoints 
was necessary. 

We focus on the case where the forcing, /, occurs at a 
wavenumber around k[ = 10. The forcing is such that the 
turbulence is subsonic and nonhclical. We consider two 
different periodic initial conditions, 



B = (cos kiz, 0, 0)Bo (nonhelical) 
and 

B = (cos kiz, sin/ciz, 0)Bq (helical), 



(12) 



(13) 



where Bo is the amplitude of the initial field. In the fully 
helical case one may expect a different decay time because 
the magnetic helicity is a conserved quantity in the limit 
of small magnetic diffusivity. For the velocity field we use 
similar initial conditions, but we do not expect this to be 
sensitive to helicity, because kinetic helicity is not con- 
served in the limit v — > 0, and would only be conserved in 
the unphysical case v = 0. 

A detailed discussion of the initial conditions may 
at first glance appear somewhat surprising, because for 
forced turbulent flows the initial conditions are normally 
forgotten after about one turnover time. This is indeed 
the case for hydrodynamic turbulence, but not for hydro- 
magnetic turbulence if the magnetic field has net magnetic 
helicity. The reason is that, regardless of the level of tur- 
bulence, the net magnetic helicity can only change on the 
resistive time scale. Our results below confirm this and 
they are indeed in agreement with earlier model predic- 
tions (cf. Blackman & Brandenburg 2002). The situation 
would be different if the initial field was bi-helical, i.e. with 
oppositely helical contributions at different scales. This 
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Fig. 1. Kinetic and magnetic energy spectra at three dif- 
ferent times for a run with Re =150 and R m = 15. 



1 http : //www.nordita. dk/data/brandenb/pencil-code 



case has been studied elsewhere (Yousef & Brandenburg 
2003). 

In Fig. ^ we show kinetic and magnetic energy spectra 
of the run with Re = 150 and R m = 15 at three dif- 
ferent times using a resolution of 512 3 meshpoints. The 
kinetic energy shows indications of a short inertial range 
in 15 < k < 40. Below the forcing scale, in 2 < k < 9, 
velocity and magnetic fields are random and (5-correlated 
in space, giving rise to a k 2 spectrum. The magnetic en- 
ergy is substantially weaker than the kinetic energy. This 
is because here the magnetic Prandtl number is small, 
Pm =0.1, and the magnetic Reynolds number is subcrit- 
ical for dynamo action. With our definition of R m the 
critical value lies around 25 (Haugen et al. 2003). The 
small scale magnetic energy is therefore maintained by 
constantly stirring the slowly decaying large scale field. 

Given that the initial large scale field depends only on 
z, it makes sense to define a mean field by averaging over 
the x and y directions. Alternatively, one might define an 
average by Fourier filtering, but this has the disadvantage 
that not all the Reynolds rules are satisfied. For example, 
the average of a product of a mean and a fluctuating quan- 
tity would not vanish. However, for all practical purposes 
our horizontal average is nearly equivalent to a projection 
onto the k = k\ Fourier mode. Indeed, the main reason for 
forcing at a large wavenumber, fcf = 10, is that we need 
some degree of scale separation. Without scale separation, 
there would be no way of distinguishing between mean and 
fluctuating fields. Since the velocity fluctuations are con- 
stantly driven via the forcing term, it would be impossible 
to measure any decay of the mean velocity. Nevertheless, 
even with scale separation there will always be a certain 
level of noise in the mean field whose energy is (fci/fcf) 2 
times smaller than energy of the fluctuations. This means 
that we can measure an exponential decay of the mean 
field only in a certain window where nonlinear effects are 
already weak, but were the noise level is not yet reached. 
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4. Results 

4.1. Decay ofu and B 

We begin by considering the decay of a helical large scale 
magnetic field and compare it with the decay of a large 
scale helical velocity field in a purely hydrodynamic simu- 
lation; see Fig. [5] Here, large scale velocity and magnetic 
fields are defined as horizontal averages over x and y; the 
result is denoted by u and B, respectively. During the time 
interval when mean velocity and magnetic field decay ex- 
ponentially, the corresponding decay rates are determined 
as 



, dlnte 2 ) 1 / 2 
A " H = dt • 



dJn(BV^ 



A B (B) = 



dt 



(14) 



In the graphs of X u (u) and Xb(B) an exponential decay 
shows up as a plateau. The magnetic field decay is ini- 
tially slow, so Xb{B) is initially not constant, but then 
it speeds up and Xb(B) reaches a plateau. The decay of 
the velocity field is immediately fast and X u (u) lies im- 
mediately on a plateau. This suggests that the turbu- 
lent magnetic diffusivity is affected by the strong initial 
field that in turn gives rise to a quenching of the tur- 
bulent magnetic diffusivity. Strong means that the mag- 
netic field strength is comparable with the equipartition 
field strength, B cq = {/iQpu 2 ) 1 / 2 . The initially strong large 
scale flow and the associated vorticity, on the other hand, 
do not and are also not expected to affect the turbulent 
viscosity and the associated decay of this large scale flow. 
For |B <C B eq , however, both u and B decay at the same 
rates, X u and As, respectively. This allows us to calculate 

v t = X u /k\ } T} t = X B jk\ 1 (15) 



where k\ is the wavenumber of the initial large scale ve- 
locity and magnetic fields. From the present simulations, 
where fcf/fci = 10, we find 



v\ « Vt = (0.8 . . . 0.9) x u ms /k { 



(forB «B e 2 q ). 



(16) 



Once |U| has decreased below a certain level (< 0.1u rms ), 
it cannot decay further and continues to fluctuate around 
0.08it rms , corresponding to the level of the rms velocity of 
the (forced!) turbulence at k = k\ (see the dashed line in 
Fig.©. _ 

The quenching of the magnetic diffusivity, rjt — rjt(B), 
can be obtained from one and the same run by simply de- 
termining the decay rate, Xb(B), at different times, corre- 
sponding to different values of B = |B|; see Fig.|3] To de- 
scribe departures from purely exponential decay we adopt 
a B-dependent rj t expression of the form 



77t(B)= %>/(! 



aB /B 2 eq ) 



(17) 



where r] t Q is the unquenched (kinematic) value of rjt, de- 
scribed approximately by Eq. I|16(). and a is a fit param- 
eter. According to Cattaneo & Vainshtein (1991) the pa- 
rameter a is expected to be of the order of the magnetic 
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Fig. 2. Decay of large scale helical velocity and magnetic 
fields (dashed and solid lines, respectively). The graph of 
u(t) has been shifted so that both u(t) and B{t) share the 
same tangent (dash-dotted line), whose slope corresponds 
to Vt = r)t = 0.86itrms/fcf • The decay of a nonhelical mag- 
netic field is shown for comparison (dotted line). 
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Fig. 3. Dependence of the turbulent diffusion coefficient 
on the magnitude of the mean field. The initial field is 
helical and corresponds to data points on the right hand 
side of the plot. R m m 20. The data are best fitted by 
a = 8 = 0.4i? m . 



Reynolds number based on the microscopic magnetic dif- 
fusivity, 



(18) 



Figure suggests that a « 0.4i? m . 

Before we discuss the effective quenching behavior of 
rjt in more detail we should note that Eq. 1)170. and in par- 
ticular the value of a, do not apply universally and depend 
on the field geometry. This is easily demonstrated by con- 
sidering a nonhelical initial field. In that case the decay be- 
comes unquenched already for B / B 2 m 1. Equation l(T7|l 
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Fig. 4. Dependence of the turbulent diffusion coefficient 
on the magnitude of the mean field. The initial field is 
nonhclical. R m m 20. The data are best fitted by a = 1, 
independent of R m . 
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Fig. 5. Dynamical quenching model with helical and non- 
helical initial fields. The quenching parameters are g = 
(solid line) and 3 (dotted line) . The graph for the nonheli- 
cal cases has been shifted in t so that one sees that the 
decay rates are asymptotically equal at late times. 



can still be used as a reasonable fit formula, but now a = 1 
produces a good fit (independent of R m ); see Fig. 0] 

In the nonhelical case there is an initial phase where 
the field increases due to the wind-up of the large scale 
field. Since we measure rjt from the decay rate of the large 
scale field, this would formally imply negative values of 
rj t . Traces of this effect can still be seen in Fig. 0] near 

2 

B /-Bgq = 1- F° r this reason our method can only give 
reliable results if \B\ ^ 0.8-B cq . In the case of a helical 
initial field, on the other hand, we have J x B = 0, i.e. 
the large scale field is force-free and interacts only weakly 
with the turbulence. In particular, there is no significant 
amplification from the initial wind-up of the large scale 
field. 

4.2. Comparison with the dynamical quenching model 

In the case of a helical field and for B / B 2 q > i?" 1 
the slow decay of B is related to the conservation of 
magnetic helicity. As discussed already by Blackman & 
Brandenburg (2002), this behavior is related to the phe- 
nomenon of selective decay (e.g. Montgomery et al. 1978) 
and can be described by the dynamical quenching model. 
This model goes back to an early paper by Kleeorin & 
Ruzmaikin (1982, see also Kleeorin et al. 1995), but it ap- 
plies even to the case where the turbulence is nonhelical 
and where there is no a effect in the usual sense. However, 
the magnetic contribution to a is still non-vanishing be- 
cause it is driven by the helicity of the large scale field. 

To demonstrate this quantitatively we solve, in the one 
mode approximation (fe = kx) with B — Bexp(ifciz), the 
mean-field induction equation 



dB 
~d7 



together with the dynamical a-quenching formula 
[Eq. (13) of Blackman & Brandenburg (2002)] 



da 
~dt = 

where 



~ Re(£ -B) 
m B 2 



£ = aB — r/t.ifci x B 



(20) 



(21) 



is the electromotive force, and i? m is defined as the ratio 
Vto/Vi which is expected to be close to the value of R m as 
defined by Eq. (JTSJ. 

In Fig. [5] we show the evolution of Bj B eq for helical 
and nonhelical initial conditions, B oc (1, i, 0) and B oc 
(1,0,0), respectively. In the case of a nonhclical field, the 
decay rate is not quenched at all, but in the helical case 
quenching sets in for B /B 2 q *b . 

2 

In the helical case, the onset of quenching at B /B 2 q m 
R^ 1 is well reproduced by the simulations. In the nonheli- 
cal case, however, some weaker form of quenching sets in 

2 

when B /B 2 q rj 1 (Fig. |3J. We refer to this as standard 
quenching (e.g. Kitchatinov et al. 1994) which is known 
to be always present. In Blackman & Brandenburg (2002) 
this was modeled by allowing in Eq. l(2"T|l 7ft to be In- 
dependent. They adopted the formula 



Th = W(l +~9\{B)\/B eq ) 



(22) 



ifei x £ - i]k\B 



(19) 



and found that, for a range of different values of i? m , g = 3 
resulted in a good description of the simulations of cyclic 
ail-type dynamos (Brandenburg et al. 2002). We empha- 
size that this rj t is not used in a diagnostic way as in 
Eq. I|17|l . but rather in the numerical solution of Eqs 119f) 
and (|20() . The resulting decay law, shown as a dotted line 
in Fig. |S1 agrees now with the decay law seen in the tur- 
bulence simulations (Fig. EJ. The helical case with g = 3 
is still compatible with the simulations. 
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Fig. 6. Decay rate for three different values of Re and 
R m = 20 (fixed), corresponding to values of P m = R m /Re 
ranging from 0.1 to 1. All three curves have a plateau 
where the value of Ab is the same. For R m = 80 and 150 
the graphs of As have been shifted in t so that all three 
graphs show the plateau in approximately the same time 
interval. 

5. Independence of microscopic viscosity 

Finally we need to show that the turbulent magnetic 
Prandtl number is indeed independent of the microscopic 
magnetic Prandtl number. In Fig. we plot the decay 
rates, obtained by differentiating In £?(£), for three differ- 
ent values of the microscopic viscosity, keeping 77 fixed. 
The resulting values of the flow Reynolds number, Re = 
u lms k[/is, vary between 20 and 150, giving P m in the range 
between 0.1 and 1. Within plot accuracy the three values 
of As turn out to be identical in the interval where the 
decay is exponential. 

The duration of this interval is u rms kfAt w 200, 
which is comparable to the time interval in Fig. [21 dur- 
ing which the decay is exponential. In one of the three 
cases (i? m = 20) the initial field was rather strong, so 
that it takes a long time before the magnetic hclicity con- 
straint becomes unimportant so that the field can decay 
exponentially (u lms kft w 800). 

The numerical resolution used in most of the models is 
128 3 mesh points. However, as Re is increased, higher res- 
olution is required. For Re=80 we used 256 3 mesh points 
and for Re=150 we used 512 3 mesh points. This implies 
mesh a Reynolds number, u rras Ax/v, based on the mesh 
spacing Ax, of about 18. Empirically we know that larger 
values are not generally possible. 

6. Conclusions 

The turbulence simulations presented here have shown 
that the turbulent magnetic Prandtl number is always of 
order unity, regardless of the values of the microscopic 
magnetic Prandtl number. Under the assumption of in- 



compressibility, both the quasilinear approximation and 
the renormalization group approach give turbulent mag- 
netic Prandtl numbers somewhat below unity, which is 
related to the pressure term in the momentum equation. 
Here we find instead P m , t ~ 1 . There are several plausible 
reasons for this discrepancy: (i) our simulations are ac- 
tually weakly compressible, (ii) they are non-steady and, 
(iii) the idealizing assumptions made in the analytic ap- 
proaches may not be justified. 

Our results have also shown that, for nonhelical mag- 
netic fields, the turbulent magnetic diffusivity is quenched 
when the magnetic energy becomes comparable to the ki- 
netic energy. For helical magnetic fields, however, an ap- 
parent suppression of the decay rate is observed which 
agrees with predictions from a dynamical quenching 
model. If this suppression is described by an algebraic 
expression, quenching would set in for magnetic energies 
much below the kinetic energy. 

The present work demonstrates that the dynamical 
quenching approach is not restricted to dynamos, but it 
can also deal with decay problems, as was already men- 
tioned in Blackman & Brandenburg (2002). The dynami- 
cal quenching model is usually formulated in terms of a, 
but for helical mean fields J and B are parallel and the 
separation into contributions from aB and % J becomes 
less meaningful. It is for this reasons that an a term ap- 
pears in the description of the decay of helical fields, rather 
than a dynamical contribution to ?7t-quenching. 

The remaining quenching of rjt that affects both he- 
lical and nonhelical fields is consistent with an algebraic 
quenching formula that is non-catastrophic, i.e. indepen- 
dent of the microscopic magnetic diffusivity. 

Although our results suggest that the turbulent mag- 
netic Prandtl number is of order unity, we cannot claim 
that it is safe to use turbulent viscosity and magnetic 
diffusivity in a simulation of the dynamo or the magne- 
torotational instability, for example, as a replacement of a 
fully resolved simulation. First of all, the functional form 
of the turbulent transport coefficients is for realistic tur- 
bulent flows more complicated and involves in practice 
tensorial rather than scalar coefficients. Numerical evi- 
dence for this has been presented elsewhere in the context 
of shear flow turbulence (Brandenburg & Sokoloff 2002). 
Furthermore, there will be additional terms such as the 
a-effect (see Sect. 14 .2|) and the AKA-effect (Frisch et al. 
1987; see also Brandenburg & Rekowski 2001). Most im- 
portantly, turbulent transport may be nonlocal, as is well 
known in meteorology when modeling atmospheric flows 
(Stull 1984, Ebert et al. 1989), where the turbulent trans- 
port is described by so-called transilient matrices (see also 
Miesch et al. 2000 for examples of astrophysical convec- 
tion). Nonlocal transport means that the transport coef- 
ficients have to be replaced by integral kernels. In Fourier 
space, the convolution with an integral kernel corresponds 
to a multiplication with a wavenumber dependent factor. 
There is indeed some evidence that the main contribution 
comes only from the smallest wavenumbers (Brandenburg 
& Sokoloff 2002). This is primarily a consequence of a lack 
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of scale separation in naturally forced turbulence, such as 
shear flows or convection. In the present context, however, 
this is not an issue because we have deliberately considered 
the case where the scale of the turbulent eddies is much 
smaller than the scale of the large scale field (kf/ki = 10). 

Finally, we wish to point out that studies of instabili- 
ties (e.g. the magnetorotational or the dynamo instability) 
using turbulent transport coefficients can sometimes lead 
to paradoxical situations. In the case of solar convection, 
for example, one expects from mixing length theory that 
turbulent viscosity and thermal diffusivity are on the order 
of a few times 10 12 cm 2 s _1 . However, using such values in 
a global model of the sun leads to an instability (Rudiger 
1989, Rudiger & Spahn 1992), which is in fact nothing but 
a repetition of the original convection instability that leads 
to turbulence in the first place (Tuominen et al. 1994). It 
is therefore plausible that the actual values of the turbu- 
lent transport coefficients should rather be close to the 
those for marginal stability. This would lead to a global 
constraint similar to the magnetic helicity constraint that 
governs the nonlinear behavior of the a-effect in helical 
hydromagnetic turbulence. At present, however, there is 
no theoretical framework that allows self-consistent mod- 
eling of convection using mean-field theory. 
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